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Abstract 

It is demonstrated that there exists a direct correlation between chemical freeze- 
out point and the softest point of the equation of state where the pressure divided 
by the energy density, p(e)/e, has a minimum. A dynamical model is given as 
an example where the passage of the softest point coincides with the condition for 
chemical freeze-out, namely an average energy per hadron ~ 1 GeV. The sensitivity 
of the result to the equation of state used is discussed. 
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Over the last few years the question of chemical equilibrium in heavy ion collisions has 
attracted much attention |jl|, ^, ^ . Assuming thermal and chemical equilibrium within a 
statistical model, it has now been shown that it is indeed possible to describe the hadronic 
abundances produced at beam energies ranging from 1 to 200 AGeV. The observation 
was made that the parameters of the chemical freeze-out curve obtained at CERN/SPS, 
BNL/AGS and GSI/SIS all lie on a unique freeze-out curve in the temperature T - 
chemical potential fiB plane 0. Recently, a surprisingly simple interpretation of this 
curve has been proposed: the hadronic composition at the final state is determined solely 
by the average energy per hadron, {Ehad) / {Nhad), of approximately 1 GeV in the rest 
frame of the system under consideration P, 0. Using this finding {Ehad) / {Nhad) ~ 1 GeV 
as a heuristic definition for the chemical freeze-out point, in this letter we study physical 
conditions which are realized around this point in evolution of the system. We show that 
the chemical freeze-out point is intimately related to the softest point of equation of state 
defined by the minimum of the pressure-to-energy density ratio, p{e)/e, as a function of 

Our considerations are essentially based on the recently proposed Mixed Phase (MP) 



model [P, 10 1 which is consistent with the available QCD lattice data |rT|. The underlying 



assumption of the MP model is that unbound quarks and gluons may coexist with hadrons 
forming an homogeneous quark/gluon-hadron phase 10|. Since the mean distance 
between hadrons and quarks/gluons in this mixed phase may be of the same order as 
that between hadrons, their interaction with unbound quarks/gluons plays an important 
role defining the order of the phase transition. Recently the importance of quark/gluon- 
hadron interactions was also discussed in the context of dilepton production . 

Within the MP model ||^, the effective Hamiltonian is written in the quasiparticle 
approximation with the density- dependent mean-field interaction. Under quite general 
requirements of confinement for color charges, the mean-field potential of quarks and 
gluons is approximated by the following form: 

UM = u,{p) = ^ (1) 

with the total density of quarks and gluons 

P = Pq + Pg + Y. ^JPj 
j 

where Pg and Pg are the densities of unbound quarks and gluons outside hadrons, while 
Pj is the density and nj is the number of valence quarks inside the hadron of type j. The 
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presence of the total density p in (HD describes the interaction between all components of 
the mixed phase. The approximation (|l]) recovers two important limiting cases of the QCD 
interaction. Namely, if p — the interaction potential goes to infinity, i.e. an infinite 
energy is required to create an isolated quark or gluon which ensures the confinement of 
color objects. In the other extreme of high energy density corresponding to p ^ oo, we 
obtain the asymptotic freedom regime. 

The use of the density-dependent potential (|lj) for quarks and the hadronic potential 
described by a modified non-linear mean-field model |jl3| requires certain constraints to 
satisfy thermodynamic consistency |P, |10[- For the chosen form of the Hamiltonian these 
conditions show that Ug{p) and Uq{p) should be independent of the temperature. From 
these requirements one also obtains an expression for the quark-hadron potential 0. In 
the case when the quark-gluon component is neglected, the MP model is reduced to the 



interacting hadron gas model with a non-linear mean-field interaction |jT3|. 

A detailed study of the pure gluonic medium with SU (3) color symmetry and a first 
order phase transition allows to fix the values of 7 = 0.62 and = 250 MeV . These 

values are then generalized to the S't/(3) system with dynamical quarks. For two light 
flavors at zero baryon density, = 0, the MP model is consistent with the results from 
lattice QCD |Tl]] with the deconfinement temperature Tdec = 153 MeV and the crossover 
type of phase transition. The model is then proposed to be extended to baryon-rich 
systems in a parameter-free way ||^. 

A particular consequence of the MP model is that for ub = the softest point 
of the equation of state is located at comparatively low values of the energy density: 
Esp ~ 0.45 GeV/ fm? which is in a good agreement with recent lattice estimates [ pJ| . 
This value of Ssp is close to the energy density inside a nucleon. Thus, reaching this value 
signals that we are dealing with a single 'big' hadron consisting of deconfined matter. 
For baryonic matter the softest point is gradually washed out at ns > 0.4 tiq. As shown 
in this thermodynamic behavior differs drastically from both the interacting hadron 

gas model which has no softest point and the two-phase approach, based on the bag 
model, having by construction a first order phase transition and the softest point at 
Ssp > ^ GeV/ frn? independent of 0. These differences should manifest themselves 
in the expansion dynamics. 

In Fig.l we show the evolution trajectories for central Au + Au collisions in the T — fip 
plane together with the freeze-out parameters obtained from hadronic abundances. Using 
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the quark-gluon string model |T^, the initial state was estimated as a state of hot and 
dense nuclear matter inside a cylinder in the center-of-mass frame with radius R = 4 fm 
and length L = 2R/'yc.m.- The calculated temporal behavior of energy and baryon densi- 
ties is close to that in RQMD or UrQMD transport models and the procedure of defining 



the initial state of a fireball is described in |T0[. The subsequent isentropic expansion 
of this fireball is treated within the Frankfurt expansion model [|16| based on a scaled 
hydrodynamical prescription assuming t ^ V. As seen in Fig.l, the turning points of 
these trajectories, i.e. the points where dT/dfiB changes the sign, are in a good agree- 
ment with the extracted chemical freeze-out parameters and clearly correlate with the 
smooth curve for the ideal gas model with {Ehad) / {Nhad) = 1 GeV defined above as a 
condition of chemical freeze-out. Existence of the turning point is related to two limiting 
equilibrium regimes for temperature behavior of the chemical potential: ultrarelativistic 
regime dominated by light particles (quarks, pions) when fiB ~ T"^ and non-relativistic 
one defined mainly by nucleons and deltas with /ib ~ T^^ 



Footnote: Other approach have been recently applied in |]14| where conditions for ther- 
modynamic equilibrium in heavy ion collisions were studied within the transport UrQMD 
model and the calculation results at a fixed time moment were approximated by the ideal 
hadron gas equation of state to extract the time dependence of T and fis- It is of interest 
that our results for dynamical trajectories below the turning point practically coincide 
with that from for equilibrium part of T — /i^ trajectory. However, at higher temper- 
ature, where according to []T1| the system is in a non- equilibrium state, their results differ 
from ours and exhibit no turning point. 



. Uncertainty in the initial e and in 10-20% shifts slightly the dynamical trajectory 
but the turning point nevertheless stays on the {Ehad) / {Nhad) = 1 GeV curve. Notice 
that there is a finite region in T where dfis/dT ^ 0, i.e. the chemical potential in the 
expanding system is kept constant around this turning point. 

The observed correlation is further elucidated in Fig. 2. The quantity p/e is closely 
related to the square of the velocity of sound and characterizes the expansion speed 



Footnote: In simplified hydrodynamic models was shown that, for example, the transverse 
expansion of a cylindrical source is governed by the pressure-to-enthalpy ratio, p/ {p + 
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e) |T^, |T8|, i.e. by p/e rather than by the sound velocity. 



, so the system hves for an appreciable amount of time around the softest point which 
facilitates to reach the chemical equilibrium. This statement is evidenced directly in the 
right-hand side of Fig.2 where the time evolution of p/e is shown: The system spends 
about 1/3 of the total expansion time in a state near the softest point. During this time 
the baryon chemical potential fiB remains practically constant. 

In Fig.2 (left-hand side) the position of the softest point correlates with the average 
energy per hadron being about 1 GeV for all beam energies except for 2 AGeV case. One 
should note that the quantity e/phad, where e is the total energy density, coincides with 
{Ehad) I {Nhad) considered in only in the case when there are no unbound quarks/gluons 
in the system. In the MP model, all components are interacting with each other and 
therefore the quantity {E^ad) is not well defined. However the admixture of unbound 
quarks at the softest point Ssp amounts to about 13% and 8% at beam energies Eiab = 150 
and 10 AGeV , correspondingly. Thus thermodynamical quantities and in particular the 
ratios of hadron abundance in the MP and resonance gas model are very close to each 
other at the freeze-out point. 

The MP equation of state plays a decisive role for the regularity considered here, 
describing both the order of the phase transition and the deconfinement temperature. 
The two-phase (bag) model exhibits a first order phase transition with T^ec = 160 MeV 
and has a spatially separated Gibbs mixed phase but the corresponding trajectories in 
the T — fiB plane are quite different from those in the MP model as shown in |T^. The 



exit point from the Gibbs mixed phase at Eiab = 150 AGeV is close to the corresponding 
freeze-out point in Fig.l. However it does not lead to the observed correlation with the 
softest point position in the whole energy range considered. The interacting hadron gas 
model has no softest point effect as was demonstrated in |T^ but nevertheless it exhibits 



the turning point in T — fiB plane. These facts are seen also from Fig.l and Fig.2 where for 
Eiab = 2 AGeV quark-hadron interactions are negligible (~ 1%). Indeed, at this energy 
instead of a minimum one has a monotonic fall-off specific for hadronic models with only 
a small irregularity in p/e near the point e/phad ~ 1 GeV 



Footnote: Note that at the SIS energies the chemical freeze-out point practically coincides 
with the thermal freeze-out [p!9[| . 
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. At higher energies the djTiamical trajectories are very different for hadronic and mixed 
phase models but the positions of turning points in the T — jis plane turn out to be still 
rather close to each other and both lying on the curve corresponding to e/ phad = 1 GeV . 
Therefore, one could expect similar results for the ratios of hadronic abundance in these 
two models. However the model predictions differ essentially as long as evolution of the 
system is concerned. 

It is noteworthy that similarly to the results presented in Fig.2, the softest point 
of the equation of state correlates with an average energy per quark, e/p ~ 350 MeV 
which is close to the constituent quark mass. So, at higher values of e/p we are deal- 
ing with a strongly-interacting mixture of highly-excited hadrons and unbound massive 
quarks/gluons forming (in accordance with Landau's idea pO|) an 'amorphous' fluid suit- 
able for hydrodynamic treatment. Below the softest point the interaction deceases, the rel- 
ative fraction of unbound quarks/gluons is getting negligible, higher hadronic resonances 
decay into baryons and light mesons. At this stage the application of hydrodynamics 
becomes questionable. 

In summary, we have shown in this paper an example of an equation of state where 
there exists a direct correlation between chemical freeze-out point and the softest point 
of the equation of state. This description correlates with the observation that chemical 
freeze-out occurs when the average energy per hadron, {Ehad) I {Nhad), drops below the 
value of 1 GeV as found in 0. In the examples considered such a clear correlation is 
observed for the equation of state of the mixed phase model, where the energy per hadron 
seems to be a relevant observable defining chemical freeze-out conditions. 
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Figure 1: The compiled chemical freeze-out parameters (from ||^, 0|) obtained from the 
observed hadronic abundances and dynamical trajectories calculated for central Au + Au 
collisions at different beam energies Eiab with the mixed phase equation of state. The 
smooth dashed curve is calculated in the ideal hadron gas model for {Ehad) / {Nhad) = 
IGeV i. 
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Figure 2: The ratio of pressure to energy density, p/e, versus the average energy per 
hadron, e/phad, (left-hand side) and its time evolution for different systems (right-hand 
side). The curve for pp collisions at = 40 GeV is calculated for isentropic expansion 
of a sphere with R = 1 fm. Other cases are calculated for central Au + Au collisions at 
the given beam energy and under the same conditions as in Fig.l. 
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